arXiv:hep-lat/9808026v2 8 Sep 1998 


KUNS-1526 
HE(TH) 98/12 
^ep-lat/9808026 

Axial vector current of exact chiral symmetry 
on the lattice 

Yoshio Kikukawa* 

Department of Physics, Kyoto University 
Kyoto 606-01, Japan 

and 

Atsushi Yamada^ 

Department of Physics, University of Tokyo 
Tokyo 113, Japan 

August, 1998 


Abstract 

We discuss the exact chiral symmetry and its spontaneous break¬ 
down in lattice QCD with the Dirac operators satisfying the Ginsparg- 
Wilson relation. The axial vector current, which turns out to be related 
to the vector current simply by the insertion of the operator 75(1 —all), 
is explicitly constructed in the case of the Neuberger-Dirac operator. 
We also consider an Euclidean proof of the Nambu-Goldstone theorem 
using the Ward-Takahashi identity of this symmetry. 
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Recently Neuberger has proposed a Dirac operator which describes ex¬ 
actly massless fermions on a lattice [Q, ^ based on the overlap formalism of 
chiral determinant il- Its explicit form is known as 

a-D„™= (1) 

where H is the hermitian Wilson-Dirac operator defined by 

^ = 75 , (0<mo<2). (2) 

Here SU{N) gauge group is assumed and is the gauge covariant forward 
difference operator. 

This Dirac operator satisfies the Ginsparg-Wilson relation i I 10,0 
Dnm'y^ ^ ^ ■ ( 3 ) 

I 


This relation guarantees that the effects of the chiral symmetry breaking 
terms in the Dirac operator appear only in local terms and physically irrel¬ 
evant. This is the clue to escape the Nielsen-Ninomiya theorem |TJ]. 

Subsequently, Liischer pointed out that the Ginsparg-Wilson relation 
implies an exact symmetry of the fermion action ||T^. The chiral transfor¬ 
mation proposed is given as 




1-^D 


75- 


( 4 ) 


This chiral transformation can be extended to the flavor non-singlet chi¬ 
ral transformations. Liischer also observed that for the flavor-singlet chi¬ 
ral transformation the functional integral measure is not invariant and its 
anomalous variant is indeed given in terms of the index of the Dirac operator. 

Locality of the Neuberger-Dirac operator in the presence of gauge field is 
not quite obvious, although it provides a completely satisfactory local solu¬ 
tion of the Ginsparg-Wilson relation in free theory |^, 151. This question has 


been examined analytically and numerically by several authors. Hernandes, 
Jansen and Liischer have given a proof of the locality for a certain set of 
bounded small gauge fields and also for the case with an isolated zero mode 

^Another Dirac operator which satisfies the Ginsparg-Wilson relation has been pro¬ 
posed by Hasenfratz et. al. i I i 0 0 based on the renormalization group method 
10,10 . As to this Hxed point Dirac operator, the authors refer the reader to the above 
cited references. 
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of H 1^. Locality in dynamically generated gauge fields has been exam¬ 
ined by same authors for SU{3) and by Montvay for SU{2) Hi- According 
to the former investigation, the Neuberger-Dirac operator seems local even 
with the gauge fields for (3 > 6.0. 

The question of the practical use of the Neuberger-Dirac operators also 
does not seem to be serious, although the Dirac operator involves the inverse 
square root of and appears at first sight to require the full storage of the 
size of the Dirac operator. In this respect, an approximate formula has been 


given by Neuberger |]T], for which the conjugate gradient method for the 
(hermitian) Wilson-Dirac operator is applicable with the help of the shift 
technique for different fermion masses ||l^ . Other approaches have also been 
investigated by Chiu and by Edwards, Heller and Narayanan |21|. 

The weak coupling perturbation theory of the massless QCD with the 
Neuberger-Dirac operator has been discussed by the authors |^2[ . The axial 
17(1) anomaly is evaluated explicitly and is shown to have the correct form of 
the topological charge density for perturbative backgrounds, supplementing 
the original calculation in |^]. 

In this paper, we discuss the axial Ward-Takahashi identity associated 
with the exact chiral symmetry in some detail. We first consider the axial 
vector current which, following the Noether’s procedure, is dehned through 
the variation of the action by the “local” chiral transformation. 


S5SF = 


nml ^ Vi Jim \ ^ J nl J 

— ® ^ ( Ipn {UnlOll OtnDnl\ T5 (f 

nml 

= XI I X i’nKff, (n, m) V’m [ • 


(5) 


Here we have used the Ginsparg-Wilson relation to derive the second line. 
As we see from this structure of the variation of the action, the axial vector 
current may be constructed following the general procedure discussed by 
Ginsparg and Wilson ||5|, leading to the third line, where is the forward 
difference operator on the lattice. ^ In a practical point of view, however, 

^Another option about the choice of the axial vector current is to use the current 
induced from the symmetric part of the Dirac operator i [D, 75 ] 75 . In this choice, as 
argued by P. Hasenfratz||li|], the explicit breaking term of the action do not induce the 
additive correction to quark mass, the renormalization of axial and vector currents and the 
mixing between four-fermion operators with wrong chirality, except the axial anomaly fl^. 
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it would be useful to have the explicit form of the axial vector current, ^ 
at least, in the case with the Neuberger-Dirac operator. We will derive the 
explicit formula of the vector and axial vector currents in this case and also 
in the case with the Dirac operator in the rational approximation of Ref. 
[18|. Then, we discuss an Euclidean proof of the Nambu-Goldstone theorem 
of the spontaneous breakdown of the exact chiral symmetry, using the axial 
Ward-Takahashi identity. 

We first note from Eq. (^) that the axial vector current has a simple 
relation to the vector current which is defined by 

SSp — CL 'y ^ 'tpn ^nDnl\ '^m 

nml 


= X! i X! m) V’n 


( 6 ) 


Once the vector current is constructed, the kernel of the axial vector current 
is obtained simply by the insertion of the operator 75(1 — aD)\ 

Kl = Ki^-l^{l-aD). (7) 

The above relation can naturally be understood if we recall the fact that 
the operator 75 (1 — aD) satisfies the identity 

{75 (1 - aD)Y = 1 , ( 8 ) 


which follows from the Ginsparg-Wilson relation, as discussed by Narayanan 
[p5| and its eigenvalues may be regarded to define another chirality analo¬ 
gous to 75, as considered by Hasenfratz, Niedermayer and Liischer [^]. This 
chirality depends on the gauge fields. In fact, if we define the chiral compo¬ 
nents of 'tp and p) through 75 (1 — aD) and 75, respectively, 

75(l-aD)V^± = iV’i, (9) 

'^L,Rlb = ( 10 ) 

the fermion action is decomposed into chiral parts as follows: 

{pDip = iPrPlDiPp + tPlPrDiP-. ( 11 ) 


This leads to the factorization of the partition function into chiral determi¬ 
nants in the overlap formalism |^]. In the kernel of the axial vector current 


®The covariant current in the overlap formalism has been discussed by Neuberger in 
and 1^. The relation between the covariant current and the Noether current obtained 
here will be discussed in a separate paper 
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of Eq. the operator 75(1 — aD) is to assign charges to V'i opposite each 
other, contrary to the vector charges. 

In this respect, we also notice the following fact. For a chiral trans¬ 
formation of the type given by Liischer to lead the exact symmetry of the 
action by the Ginsparg-Wilson relation, the weights of the Dirac operator 
in the transformation for ijj and ip should sum up to unity, but otherwise 
may be arbitrary. We can easily see by considering “local” transformations 
following the Noether’s procedure, that all such chiral transformations lead 
to the same axial vector current as given by Eq. i) and the associated ax¬ 
ial Ward-Takahashi identities are physically equivalent. For example, the 
following chiral transformation also leaves the action invariant 

= Xl'lS (1 - <5^n = ^n75, (12) 

m 

and gives more natural definition of chiral transformation in view of the 
chiral decomposition Eqs. and (11^), as discussed in [^. We will con¬ 
sider this version of the exact chiral symmetry later in our discussion of the 
spontaneous chiral symmetry breaking. 

This particular chiral symmetry, in the case of the Neuberger-Dirac op¬ 
erator, can be regarded as a direct consequence of the symmetry of quantum 
Hamiltonian in the underlying overlap formalism |^]: 

'hio — ® ^ ® ^ ( 1 ^) 

nm nm 

Wd5 = (14) 

n n 


This can be seen from the relation between the Dirac field variables and 
-ipn and the grassman representatives of the canonical variables (a„,,al) and 

{bn,bi) 113: 



Then the symmetry transformation of the canonical variables 

6an = -an, 5ai= ai, 6bn = bn, 6bl = -b'l, (17) 


in fact induces Eq. (1^). 
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Now let us construct the vector current in the case of the Neuberger- 
Dirac operator. The covariant transformation of the Dirac operator under 
the local U{1) transformation 

Dnm Dnm (18) 

is induced by the local transformation of the overlap Hamiltonian H 

U _> — p-iotn TT Jam ('I O'! 


as 




H' 


m 


( 20 ) 


Then, at the first order in the infinitesimal transformation parameter a, we 
have 


O'Dnm Olm Oin — ^5 

= 75 


{(iit), 






SH={Hkiai-akHi,i} 


where 


SH = Hnmam-anHnm = '^d^aiWi^{n,m), 


( 21 ) 


( 22 ) 


^) 2 ^nlSl+ji,m Ulfj, + 2 {'y^l T 1) ^n,l+ji5lm • 


(23) 


Using the integral representation, 


we obtain 



I _ dt 1 

^/}p 7-00 TT + i72 ’ 


(24) 


5H 


1 


y/W 

°° dt 


-H 

1 


TT 


t 2 + i72 


{tHH - HSHh) 


1 

WTiP 


(25) 
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Therefore, the kernel of the vector current is given by 


aKi^{n,m) = 75 |y 


dt 


-oo vr 


- HWi„H 




(26) 

It also follows from the Ginsparg-Wilson relation that the kernel of the 
vector current satisfies an identity given by 


Kfj. • 75 (1 - aD) + (1 - aD) 75 • = 0. (27) 

Accordingly, using the relation Eq. (0) , the axial vector current is ob¬ 
tained as 

Am = (28) 



The flavor non-singlet currents can be obtained by including the generators 
T“ of the rotation in the flavor space in the above expression. 

Here we comment on the factorization property of the vector and axial 
currents into left- and right-handed currents. The vector and axial vector 
currents in Eqs. (I) and (^) are defined by introducing the local transfor¬ 
mations 


d'lpn — 

for vector current and 


(30) 


d'tpn = 75(1 - aD)nmi^m, d'tpn = anV’n75, (31) 

m 

for axial vector current, respectively, and the kernels of these currents thus 
constructed respect the simple relation Eq. (Here we have adopted the 
chiral transformation Eq. (|T^ ) rather than Eq. (Q) for the reason explained 
before.) In order to obtain the right-handed chiral current, for example, we 
may combine these two transformation as 


where 


6ll^n — Ctn i^+)nm 

m 

— ^n'4^nPLi 

(32) 

„ 1 ± 75 ( 1 - 0 - 0 ) 

2 

P 1 ±75 

Pr,l - 2 • 

(33) 
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Then we obtain the right-handed current as 


~ 9 ~ X] ■ P+)nm 


( 34 ) 


In the point of view of the chiral decomposition of Eqs. (I), (p^ and ([Il|), 
however, this current is not completely chiral since it contains not only the 
right-handed components but also This is because P± do not 

commute with the local variation operation (or operator 5 ) and the above 
local transformation is not restricted to the right-handed chiral component 
{ipR, 'ip+) only. To obtain the factorized chiral current which consists of only 
the right-handed components (V’+,'(A_r), we may consider the right-handed 
local chiral transformation given as follows: 


S'tpn = (P+)nl + '^rn, 5lpn = -an'lpnPl- 


(35) 


ml 


Here P+ in front of the local parameter is necessary so that the change 
of the variable along the chiral transformation is consistent with the chiral 
decomposition: 

1pn+ = i’n+ + S'lpn- (36) 

Then we obtain the factorized right-handed current 




'tpn 


(37) 


The transformation Eq. (35) gives rise to half of the covariant anomaly as 
the Jacobian factor of the path integral measure. The local transformation 
for the other component (V’-,'0 l) and its current are obtained in the same 
way. As it should be, the axial vector like combination of these local trans¬ 
formations is reduced to Eq. ( 0 ) when the transformation parameter is 
space independent. These currents are easily generalized to the case of the 
flavor non-singlet chiral transformations SU{Nf)R x SU{Nf)L. 

As we can see from Eq. (15), in the case of the Neuberger-Dirac opera¬ 
tor, the local transformation Eq. (EsI) can be regarded to be induced from 


the local transformation of the canonical variable (a^, a)j) in the overlap for¬ 


malism. Therefore the covariant and chiral (right-handed) current Eq. (^) 
corresponds to the covariant current in the overlap formalism [23, 271. 

Our expression of the axial vector current given by Eqs. (|^), ( p8| ) and 
( ^^ contains an integration over the parameter t. Eor analytical studies, 
for example, the weak coupling expansion, this integration can be explicitly 


performed j^j. In numerical applications, however, we would need to treat 
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it in some approximation: the integration may be replaced by a certain 
discrete summation. Such an approach is closely related to the rational 
approximation of the Dirac operator considered by Neuberger |]^.^ In fact, 
if we change the variable t to the variable which has a compact region as 


TT 

t = tan —cr, a G [—1, IJ, 

the integral representation Eq. (p4|), for example, may be rewritten as 

1 _ dt 1 _ D 1 1 

y/Jp 7-00 vr Jo ^cos"^ |(T tan^ |cr + 11“^ 

By discretizing the variable cr as 



(38) 


(39) 


(40) 


we are lead to the rational approximation to the Neuberger-Dirac operator. 




1 + 75 



(41) 


and 



(n) 



Ms-\) tan2 - \) + 


(42) 


Similarly, we may replace the integration over the parameter t in the expres¬ 
sions of the current kernels to the discrete summation. The same results are 
obtained using the expressions Eqs. ( ^T|) and (|4^ ) from the beginning. Thus 
we obtain the axial vector current in the case with the Neuberger-Dirac 
operator in the rational approximation as follows: 


aK, 


in) 

Ifl 


K] 


//i 


= 


^ ^{s-\) tan2 ^{s - \) + 


tan2^(s- 1) +dd2- 


5(n) 


= K, 


(n) 

1{1 


H 




ViPJ 


(43) 

(44) 


^We would like to thank H. Neuberger for suggesting this relation. 
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It may also be possible to use the expansion in Legendre polynomials for 
the inverse sqnare root of as adopted by Hernandes, Jansen and Liischer 
[p8| , [l6| . With the help of the integral representation Eq. we can derive 
the expansion of the current kernels in terms of Legendre polynomials. 

Finally, we consider the axial Ward-Takahashi identity associated with 


the flavor non-singlet chiral transformation Eq. (1^), 


Hn = X! ^“75 (1 - Xljm, S-lpn = V’nTsT'" 


(45) 


In particular, we discuss the chiral multiplet of the scalar and pseudo-scalar 
density operators considered by Neuberger, Hasenfratz and Chandrasekha- 


ran 


S (V’mTsV’m) = 2 

7 \ ^ / ml 


(46) 


Its axial Ward-Takahashi identity is given as follows; 


a; + 6nm6^^ E (^ " f ^ (^7) 


Here we assume that the vacuum expectation values are defined by the 
functional integrals over both the fermion and the gauge fields. In the case of 
the Neuberger-Dirac operator, we also assnme that the gauge fields satisfies 
the smoothness condition discussed in |T^] by the choice of an appropriate 
gauge held action. Note that, with this constraint, the kernel of the axial 
vector current can be shown to be exponentially (uniformly) bounded [p^]. 

Then the identity implies that the massless pole would appear in the 
correlation function of the axial vector current and the pseudo scalar density 

F.T. (p) = 

(48) 

if the spontaneous breakdown of the exact chiral symmetry occurs 

(49) 
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This consists an Euclidean proof of the Nambu-Goldstone theorem on the 
lattice. Let us assume the coupling of the Nambu-Goldstone boson to the 
axial vector current and the pseudo scalar density as follows: 

(0|^nT“75V'„|vr'(p)) = 

Then we obtain the relation, 

F^zl = 

which is well-known in the continuum theory. 

Comments about the choice of the multiplet of the scalar and pseudo 
scalar density operators are in order. Under the “global” chiral transforma¬ 
tion Eq. ( 0 ). it is easily verified that the following scalar and pseudo scalar 
density operators consist an exact chiral multiplet 

Sn = 1pn'^(l-^D] Pn = - ^F)\ V'm, (53) 

77 X \ ^ / nm jYi V ^ J nm 

as 

6Sn = 2Pn, 6Pn = 23^- (54) 

Then one might think of the axial Ward-Takahashi identity which involves 
the pseudo scalar density Pn- The above chiral property, however, is not 
reflected directly in the axial Ward-Takahashi identity, because it is derived 
through the “local” chiral transformation. In fact, what we obtain in this 
case is as follows: 


(50) 

(51) 

(52) 


f)* IA P 

^fj, ^ n 


)+5 nm {Sm) + { 


l-^D 


i 


75 


1 - -D 
2 


nl 


il^i) = 0. 


(55) 

The third term in the l.h.s. is not proportional to 6nm times the chiral con¬ 
densate, contrary to Eq. (U). As far as the axial Ward-Takahashi identity 
for the multiplet of the scalar and pseudo scalar densities is concern, the 
choice of Chandrasekharan gives rise to the simple identity. 

In summary, we have considered the exact chiral symmetry and its spon¬ 
taneous breakdown in lattice QCD with the Dirac operators satisfying the 
Ginsparg-Wilson relation. We have constructed explicitly the axial vec¬ 
tor current in the cases with the Neuberger-Dirac operator and its rational 
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approximation, which turns out to be related to the vector current sim¬ 
ply by the insertion of the operator 75(1 — aD). Then we have considered 
an Euclidean proof of the Nambu-Goldstone theorem with the axial Ward- 
Takahashi identity associated with the “local” chiral transformation. The 
relation given by Eq. (^) or the explicit formula of the axial vector current 
may be useful in the numerical evaluation of the pion decay constant TV. 
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